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^ , Abstract. Possibility to detect DCC fluctuations is discussed. It is shown that 

^1 ' interactions with quark background and dissipative effects due to interactions 

^vq . in the chiral field may result in damping of fluctuations. Since the magnitude 

^SJ ' of fluctuations depends strongly on the initial state and speed of chiral phase 

^^ , transition accurate evaluation of all modifying processes is required to predict 

r^ ' observability of DCCs. 

(~| , The possiblity of producing quark-gluon plasma in relativistic heavy ion col- 

lisons is exciting especially from the point of view of observing the chiral and de- 
confinement phase transitions. Rajagopal and Wilczek have suggested that if chiral 
J^ ' restoration is of second order, nonequilibrium dynamics can generate transient do- 

J^ I mains in which macroscopic pion fields develop. As the chiral field relaxes to the 

true vaccum in such a domain, it may lead to coherent emission of pions [ |l|. This 
kind of phenomenon is called a disoriented chiral condensate (DCC). Bjorken and 
others pointed out that DCC can lead to fluctuations in the charged and neutral 
pion spectra [ Q] . If wc can observe this phenomenon we can say the quark-gluon 
plasma has been produced, and we see the signs of its hadronization. 

The possibility to resolve fluctuations of neutral and charged pions critically 
depends on the size and energy content of these domains. If the domains are roughly 
pion sized, the effect of DCC is too small to be resolved in experiments. If domains 
are large enough, the effects of DCC can have measurable consequences. All these 
signals: the isospin fluctuations, the enhancement of the pion spectrum at low px, 
and the suppression of HBT correlations are characteristics of any large coherent 
source [ |j. 

Most dynamical studies [ 0, @, @, M, 0] of DCC have been carried out within 
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the framework of the linear sigma model. Using this model, one can study the 
possibility of DCC formation and the domain size of DCC. 

The Lagrangian density of the linear sigma model can be written as : 

£ = £, + CsM =^a + \ i9f.<^ 9^<^ + d^if d^n) - U{cj, tt), (1) 

where 

U{a,^) = \{a' + n^ -v^f -Ha+(rnl!l-^y (2) 

is the so-called Mexican Hat potential, and Cq is the contribution of quarks, 

i^q = q [iltJ.d^ - 9(0- + i-fsTTf)] q. (3) 

Here q stands for the light quark fields, {u,d), while a and 7? — (7ri,7r2, tts) are 
the scalar meson and vector meson fields which together form a 4-component chiral 
field $ = {a, w). The last constant term is included to secure that the minimum of 
the potential is Umin = 0. In the linear sigma model, the fields a and 7? are treated 
independently without any constraint, and coupled to the fermion field as in the 
above equation. If we impose the constraint cr^ +7?^ = /^j '^6 will have the nonliner 
sigma model. Without the term Ha this Lagrangian is invariant with respect to 
the 5'C/l(2) (g) SUb.{2) chiral transformations. 

The parameters in this Lagrangian are chosen in such a way [ 0| that in normal 
vacuum at T = chiral symmetry is spontaneously broken and expectation values 
of the meson fields are 

(a) = /., (7r)=0, (4) 

where /7r=93 MeV is the pion decay constant. In the vacuum with broken chiral 
symmetry pions represent soft "azimuthal" excitations of the chiral field. To have 
the correct pion mass in vacuum, m^=138 MeV, one should take 



V 



2 _ (-2 '"-TT u _ f „,2 



/. - ^, H^ UK. (5) 



The parameter A is related to the sigma mass, m?^ — 2Xf^+ni^, which can be chosen 
to be about 0.6 GeV (then A « 20). Sigmas represent stiff, "radial" excitations of the 
chiral field. The remaining coupling constant g can be fixed by the requirement that 
the effective quark mass, (m^ = g{a^ + tt^)), in broken vacuum, rUq = g{a) = gU, 
coincides with the constituent quark mass in hadrons, about 1/3 of the nucleon 
mass m = mj\[/3. This gives g ~ {mM/i)/ f-,^ ~ 3.3. 

When we take the mean-field approximation, ignore all loop contributions and 
consider a and tt as classical fields, the equation of motion of chiral field can be 
written as : 

5,,5^cr(a;) -I- A [a'^{x) + n'^{x) - v"^] a{x) -H = -gps{x), 
df,d''7f{x) + X[cr^{x)+TT^{x)-v^]Tr{x) = -gpp{x). 

(6) 
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Here ps — (qq) and pp = i{q'y^fq) are scalar and pseudoscalar quark densities, 
which should be determined self-consistently from the motion of q and q in back- 
ground meson fields. The scalar and pseudoscalar densities can be represented as 

Ps{x) = ga{x)cj{x), pp{x) = ga{x)Tr{x), (7) 

where a{x) is expressed in terms of the momentum distribution of quarks and an- 
tiquarks f{x,p), 

«(^) = ^J 7^ 2J(P>M - m'{x))f{x,p) 



Here Vq is the degeneracy factor of quarks. Let us assume that the quark distribution 
is an ideal Ferimi distribution at the freeze-out time tq (tq ~ 7 fm/c) and the 
subsequent distribution for time r > tq is given by the solution of the Vlasov 
equation [ |9[]. The integral of this source term a{x) in the equation of motion is 
evaluated for the /i = 0, niq — case analytically earlier in ref. [|9[]. For ruq ^ 0, 



2Vq 1 [^ q . \-a^ ^ 

aix) = - — ^ , / , arcsm , / = — s-ag, 

^^ (27r)2a ,/T3^ 7o eN'i'+'-'< + 1 \/l + ^ 

where a = to/t^ m = in{T) and mo = ?ti(to) are the quark masses at r and tq. 
This integral can be numerically evaluated by using nine point Laguerre integral 
formula. This way we can follow the solution of the equation of motion numerically 
from arbitrary initial condition at tq during the whole symmetry breaking process. 
In ref. [ |9[ only an estimate was given for the initial growth rate of the instability. 
Considering an A^-dimensional uniform scaling expansion, [N = 0,1,2,3) the 
differential operator, d^d^, in terms of the variables, proper time, r, and appro- 
proate N-dimensional rapidity, 77, takes the form 



T T^ ' 



(9) 



where the term, —Nr^^dr, describes the damping caused by the collective expan- 
sion of the system. A similar damping can occur due to interactions, particularly 
due to thermal interactions if our system is interacting with heat-bath. To describe 
such damping we might add an additional damping term which is present even in 
the absence of collective expansion (iV = 0): —jdr- 

Assuming further that our initial condition satisfies the scaling symmetry, i.e., 
it depends only on t but not on 77, such an 77 independence is conserved by the 
equations of motion. Thus for the studies we describe in the following we drop the 
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rapidity dependent terms in the equations of motion. The arising equations in the 
(a, 7?) coordinates are 

a> = (-^dr-jadr-X{<J^+7i^) + XvAa +H -gps, 

d^TT = f-^dr--f-,dr~X{(j^+Tr^)+XvATr -gpp. (10) 

In the above equation, 77^ and 70- are damping paramters of tt and a fields. 

In estimating thermal damping of small amplitude fluctuations we may take 
advantage of the estimated viscosity coefficient of the quark-gluon plasma [ [l^] 



which was used earlier in estimating thermal nucleation rates [ 11 . These yield a 
damping of about 



7 w 0.2 - 10 -^ r = 50 - 2000 MeV (11) 

fm 



for T w 100 MeV. 



rria 


a 


h 


MeV 






400 


0.0307 


3.706 


600 


0.0526 


3.446 


800 


0.0700 


3.547 



Table 1. Fit parameters a and b to estimate the pion-width as given in eq. ([1^), 
for three sigma masses including contributions of one- and two-loop self-energy 
diagrams, based on ref . [ ^ . 

In finite temperature field theory we can estimate the damping of a field using 
the imaginary-time propagators. The damping or width is characterized by the 
imaginary part of the self-energy. These hadronic thermal widths should be inter- 
preted as the damping coefficients of wave packets propagating through a dispersive 
medium. 

In the linear cr-model one- vertex one- loop graphs do not contribute an imaginary 
part to the self energy [ |l3| thus the first contributions are given by two- vertex one- 
loop diagrams [ . The damping of the 7r-field was evaluated in [ |2| by including 
a two-vertex one-loop self-energy graph with a dressed sigma-meson propagator, 
plus two two-loop graphs and two two-loop four-vertex graph contributions to the 
pion self-energy. 

The resulting pion width can be approximated with 

r. = {aTf (12) 

where F^ and T are measured in [MeV] and the parameters a and h are given in 
Table 1. 
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For T = 100 MeV and m„ = 600 MeV, this yields F^ « 300 MeV, which 
is comparable with the value above based on analogies of damping of waves in a 
viscous fluid. 

At the last moment before quark freeze-out we have a thermal distribution both 
for the quarks and the background fields. This moment was estimated in ref. [O] to 
be around tq = 7 f m and To « 130 MeV. As an illustration we used the same initial 
condition as ref. [ p| in our studies presented here. 

cr(ro) = , 7r(To) = 0, 

CT(ro) = IMeV/fm , 7fi(To) = 5MeV/fm. (13) 




15 25 35 45 

Fig. 1. Proper time evolution of a and tti fields for H = 0, without quark-source 
terms and without damping. 



Note that this initial condition is purely radial excitation of <I> from the origin 

and thus it leads to correlated a and tt fluctuations. 
150r 




Fig. 2. Proper time evolution of a and Tr/ fields for H ^ 0, without quark-source 
terms and without damping. 



For the above initial condition, we calaulate the proper time evolution of sigma 
and pion fields. The features of different scenarios are shown in Fig. |l], Fig. g, Fig. 
I and Fig. I 
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Fig. 3. Proper time evolution of a ar 
terms but without damping, 
150 



35 45 

1 fields for H ^ 0, with quark-source 




Fig. 4. Proper time evolution of a and tti fields for H ^ 0, with quark-source 
terms and with damping, ^^ — 5.0c/fm, 7^^ = 0.3c/fm . 



In Fig. |l|, we solve the equation of motion without symmetry breaking term, 
H, and do not consider the effect of quark source terms and damping. In this case, 
we exclude the background field. The sigma and pion fields do not couple to the 
quark field and they have complete symmetry. The oscillating frequencies of sigma 
and pion fields are exactly the same. They oscillate around different values due to 
different initial conditions and the motion is radial in the (a, tti) plane. 

In Fig. g, we include the symmetry breaking term, H. The dynamics is quite 
different compared with Fig. ||. In Fig. g the sigma field oscillates nonlinearly 
around /^, and the pion field oscillates around zero, their frequencies are also dif- 
ferent. The sigma field takes about 2 fm proper time to grow from zero to reach 
the true vacuum expectation value /tt. On the other hand the pion field eventually 
tends to zero when the proper time gets large enough. 

Compared to Fig. g, we include the quark source terms in Fig. |3[ The evolution 
of sigma field does not change too much and the oscillating amplitude of pion field 
is smaller due to the source term. 
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li H ^ the transition is gradual, but fast. The source term leads to a modifi- 
cation of the dynamics. Compared to the earlier dynamical studies with a quenched 
initial condition, however, the difference is not large and we get almost the same 
rapid rollover transition. 

In Fig. El we include the thermal damping term. The sigma field takes a 
longer time to grow from zero to the true vacuum, then quickly tends to /^r almost 
without oscillation. The fluctuation of the pion field is damped out quickly. It takes 
about 20 fm proper time for the pion field to get to zero. For the observability it is 
important that the fiuctuations are rapidly damped out. 

If we add the thermal damping term in the equation of motion, we can see that 
the decrease of the amplitude of field oscillation is very fast. It is a question if 
this later one is a realistic approach, because for small systems both damping and 
fluctuations should be considered simulataneously 
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